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Abstract — Non-minimum phase characteristic of hydro 
turbine shows an opposite initial power surge in the event 
of frequency disturbance, thus possesses widely different 
characteristic than thermal generating unit. This paper 
presents Automatic Generation Control (AGC) of an 
interconnected two-area multiple-units hydro-hydro 
power system in restructured electricity market. The step 
load perturbation to such a system results in heavy 
frequency oscillations and the system is unable to regain 
its stable state as the positive real parts of some 
eigenvalue pairs confirm the inherent dynamically 
unstable characteristic of the system. To stabilize the 
frequency oscillations, impacts of Superconducting 
Magnetic Energy Storage (SMES) placed at terminal of 
area and Thyristor Controlled Phase Shifter (TCPS) 
located in series with tie-line have been investigated. Two 
cases such as (a) SMES-SMES coordination and (b) 
TCPS-SMES coordination are evaluated to compare their 
effectiveness to suppress the frequency oscillations. The 
parameters of TCPS, SMES and integral gains of AGC 
loop are optimized through craziness-based particle 
swarm optimization algorithm in order to have the 
optimal transient response of the system under different 
PoolCo and bilateral transaction in restructured 
electricity market. 

Index Terms— AGC, Deregulation, SMES, TCPS 

I. Introduction 

Automatic generation control has been used for 
several years to meet the objective of maintaining the 
system frequency at nominal value and the net tie line 
power interchange from different areas at their 
scheduled values. The concept of conventional AGC is 
discussed in [l]-[3]. Controlling the frequency has 
always been a major subject in electrical power system 
operation and is becoming much more significant 
recently with increasing size, changing structure and 
complexity in interconnected power systems. 

In a competitive electricity market, there will be 
many market players, such as generating companies 
(Gencos), distribution companies (Discos), 
transmission company (Transco), and system operator 
(SO). For stable and secure operation of a power 
system, the SO has to provide a number of ancillary 
services. One of the ancillary services is the "frequency 



regulation" based on the concept of the load frequency 
control. A detailed discussion on load frequency 
control issues in power system operation after 
deregulation is reported in [4]-[8]. 

The AGC in a restructured electricity market should 
be designed to consider different types of possible 
transactions such as PoolCo-based transactions, 
bilateral transactions, and a combination of these two 
[7]. In this new paradigm, a Disco can contract 
individually with a GENCO for power and these 
transactions are done under the supervision of the 
System Operator (SO). Hence, situations may arise 
where DISCO make contract with only hydro 
generating units in competitive market place. In this 
paper we have analyzed the issue of AGC in 
restructured power system where contracts of DISOC 
exist only with GENCO based on hydro generating 
units. 

II. Linearized Model Of An Interconnected Two- Area 
Multi Units Hydro-Hydro Power System 

Fig. 1 shows the schematic of an interconnected 
two-area power system. Each area has two GENCOs 
based on hydro generation and two DISCOs. Let 
GENCO 1, GENC02, DISCO 1, and DISC02 be in area 
I and GENC03, GENC04, DISC03, and DISC04 be 
in area II. Fig. 2 shows a complete block diagram of an 
interconnected two-area multi-units all-hydro power 
system for AGC study. The nominal system parameters 
are given in appendix [9]. The closed loop system in 
Fig. 2 is characterized in state space form as 

x = A cl x + B cl u 

where x is the state vector and u is the vector of power 
demands. State matrix a c1 is constructed from Fig. 2 
and given in appendix. Eigenvalue analysis of the state 
matrix a c1 for a given operating conditions shows the 
positive real parts for some eigenvalue pairs (shown as 
shaded in Table 1) which indicate that the system itself 
is inherently unstable. Thus, load frequency 
stabilization can never be possible for such an all-hydro 
test system after the load . 
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Fig 1 Schematic of two-areasystem in restructured power sytem 
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Fig 2 Area frequency response of system after 10% step load 
TABLE I 
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Fig. 3 shows the transient response of the system 
after 1 0% step load perturbation in area 1 . It has been 
observed that under the occurrence of load changes, a 
system frequency is heavily perturbed from its normal 
operating point and system does not regain its stable 
state. To overcome this inherently unstable situation 
and to make the system stable against any load 
disturbance, the impact of Thyristor Controlled Phase 
Shifter (TCPS) [10] and Superconducting Magnetic 
Energy Storage (SMES) [11] are analyzed. The 
structure of TCPS and SMES as frequency stabilizers 
are shown in Figs. 4 and 5 respectively. TCPS is 
located in series with tie-line near 



area 1 whereas SMES is placed at the terminal of 
area(s). 
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Fig 4: Structure of TCPS as frequency stabilizer 
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Fig 5: Structure of SMES as frequency stabilizer 
III. MATHEMATICAL PROBLEM FORMULATION 

The objectives of AGC are to reestablish primary 
frequency regulation, restore the frequency to its 
nominal 

value as quickly as possible and minimize the tie- 
line power flow oscillations between neighboring 
control areas. In order to satisfy the above 
requirements, gains (K IP K I2 ) of integral controller in 
AGC loop, parameters of TCPS ( k tcps' T tcps) and the 
parameters of SMES(K S mes> t smes> t i> T 2> T 3 andT 4 ) ) are to 
be optimized. In the present work, an Integral Square 
Error (ISE) criterion is used to minimize the objective 
function defined as: 

FigureofDemeritfFDM^yrA^+Afl+APyAT 

(1) 
where AT = a given time interval for taking samples, 

Af = incremental change in frequency and AP tie = 

incremental change in tie-line power. The objective 
function is minimized with the help CRPSO based 
optimization techniques. 

IV. Craziness Based Particle Swarm Optimization 

The PSO was first introduced by Kennedy and 
Eberhart [12]. It is an evolutionary computational 




Fig 3 Linearized model of an interconnected two-area restructured power system 
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model, a stochastic search technique based on swarm 

intelligence. 

Velocity updating equation: 



..*+i 



+ dD/fl pBestf-jc* +c2D/E2 gBest-xf (2) 



Position updating equation: 



k+i 



:xf+vf +1 (3) 



The following modifications in velocity help to 
enhance the global search ability of PSO algorithm as 
observed in CRPSO [13]. 

(i) Velocity updating as proposed in [13] may be stated 
as in the following equation: 

v k+l = r 2Q; i +(l-r2) QEOCK (pBest.-xM 

(4) 

+(l-r2) -c> (1-rl) "(gBest-x;) V } 

Local and global searches are balanced by random 
number r2 as stated in (5). Change in the direction in 
velocity may be modeled as given in the following 
equation: 



,k+\ 



rl Dslfen(r3) 

k 



vf+fl 



+ l-rl) 



rl) M 
cl -(l-rl) 



rl 



gBest - 



(5) 



sign(/3) : 



Pr(r4) = 



sign(;-4) : 



■ p Best- 
In (5), sign (r3) may be defined as 
-1 (r3<0.05) 
1 (r3 > 0.05) 

(ii) Inclusion of craziness: Diversity in the direction of 
birds flocking or fish schooling may be handled in 
traditional PSO with a predefined craziness probability. 
The particles may be crazed in accordance with the 
following equation before updating its position. 
vf +1 = vf +1 + Pr(r4) UMgn(rA) v™ ss (6) 

where, Pr(r4) and sign(r4) are defined respectively as: 

[1 H</^;) (7) 

1 (M>0,5) (8) 

-1 (r4<0.5) 

During the simulation of PSO, certain parameters 
require proper selection, as PSO is much sensitive to 
the selection of input parameters. The best chosen 
maximum population size = 50, maximum allowed 
iteration cycles = 100, best Pcraz = 0.2 (chosen after 
several experiments), best values of cl and c2 are 1.65. 
The choice of cl, c2 are very much vulnerable for PSO 
execution. The value of v ^ raziness lies between 0.25 and 

0.35. The novelty of CRPSO lies with the fact of its 
faster convergence to the true optimal solution as 
compared to its other counter parts. 

V. Simulation Results And Discussion 

The concept of a "DISCO participation matrix" 
(DPM) is used to make the easier visualization of 
contracts between GENCOs and DISCOs [7]-[8]. Each 
area is having two GENCOs and two DISCOs. Let 
GENCOl, GENC02, DISCO 1 and DISC02 be in 



areal and GENC03, GENC04, DISC03 and DISC04 
be in area2. Unlike in the traditional AGC system, a 
DISCO asks/demands a particular GENCO or 
GENCOs for load power. Thus, as a particular set of 
GENCOs are supposed to follow the load demanded by 
a DISCO, information signals must flow from a DISCO 
to a particular GENCO specifying corresponding 

demands. The demands are specified by cpf s (elements 

of DPM) and the pu MW load of a DISCO. These 
signals which were absent in traditional AGC will carry 
information as to which GENCO has to follow a load 
demanded by which DISCO. apf t (/ = 1,2,3,4) are the 

Area Control Error (ACE) participation factors of 
different GENCOS. 

~cpf\\ c P f\2 cpfn cp/W 

cpf 21 cpf 22 cpf 21 cpf '24 
cpf 11 cpf 12 cpf 11 cpf 14 

cpf 4 1 cpf 4-2 cpf 41 cpf 44 _ 

For TCPS-SMES coordination, TCPS is considered 
in series with a tie-line near area 1 and SMES is placed 
at the terminal of area2 whereas for SMES-SMES 
coordination, SMES units are considered at the 
terminal of both areas. 

Case 1: PoolCo based transaction: 

An event is simulated in which a system shown in 
Fig. 2 is subjected to step load perturbations as given in 
(10)[7]-[8]att = 5sec. 

The GENCOs in each area participate equally in 
AGC, indicated by Area Control Error (ACE) 
participationfactors opf i ( i = 1, 2, 3, 4) . 



DPM = 



(9) 



DPM = 



0.5 0.5 0.0 0.0" 




"Aftf 




"0.05" 


0.5 0.5 0.0 0.0 




^h 




0.05 


0.0 0.0 0.5 0.5 


■ 


&L3 




0.05 


0.0 0.0 0.5 0.5 




*'}.! 




0.05 






(10) 



apfi =Q.5,apf 2 =0.5,^3 = 0.5,<p/ 4 =0.5 

The 0.05 pu load is demanded by DISCO 1 and 
DISC02 identically from their local GENCOs in areal 
only, i.e. GENCOl and GENC02, which is indicated 
by each non-zero participation factor of 0.5 in DPM. 
Similarly, DISC03 and DISC04 demand 0.05 pu load 
from their local GENCOs in area 2, i.e. GENC03 and 
GENC04. As there is no contracts between inter-area 
GENCOs-DISCOs, the inter-area participation factors 
of DPM (shown shaded in (10)) are zero. This type of 
overall transaction is called as PoolCo transaction. The 
un-contracted load demand ^z, U ncoti and ^z, U ncot2 of 
the areal and area2 respectively are zero. The transient 
response of the system for the transaction given by (10) 
is shown in Fig. 6. It is clearly observed that the 
coordination of SMES-SMES and TCPS-SMES 
effectively suppress the area frequency and tie-line 
power oscillations and stabilize the system 
effectively. .The optimized parameters for both 
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coordinated system are given in Table II. Power The mutual scheduled tie-line power flows among the 
generation response shown in Fig. 6 indicates that the areas can be represented as [8]: 
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Figure 7. Case 5 



Figure fc. Case i 



initial surge of power is in the opposite direction to that 
desired due to non-minimum phase characteristic of 
hydro turbine [9]. Hence, issue of LFC with only hydro 
generating unit has become very critical and needs a 
special coordinated control of TCPS-SMES or SMES- 
SMES. Since there are no contracts of power between a 
GENCO in one area to a DISCO in another area, the 
scheduled steady state power flow over the tie line is 
zero. 

Case 2: Bilateral Contracts: 

The nonzero participation factor of DPM of (11) 
represents both PoolCo and bilateral types of 
transactions unlike that indicated by DPM of (10). The 
apfi of GENCO 1 is the highest; the apf 2 of GENC02 is 
the lowest which indicates that the participations of 
GENCO 1 and GENC02 in zeroing ACE are the 
highest and the lowest respectively. The total 
generation required of individual GENCOs can be 
calculated as [8]: 
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Table II 

Optimized Parameters obtained by CRPSO for TCPS- 

SMEs and SMES-SMES Coordination 
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The dynamic responses for combined PoolCo and 
bilateral transactions are shown in Fig. 7. In this case 
also, the coordination of SMES-SMES and TCPS- 
SMES works satisfactorily. As the inter-area bilateral 
contacts exist between the DISCOs and GENCOs, 
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scheduled tie-line power are observed on the tie-line in 
steady state also. 

Case 3: Contract violation 

It may happen that a DISCO violates a contract by 
demanding more power than that specified in the 
contract. This excess power is not contracted out to any 
GENCO. This un-contracted power must be supplied 
by the GENCOs in the same area as the DISCO. It 
must be reflected as a local load of the area but not as 
the contract demand. Consider case 2 again with a 
modification that DISCO demands 0.05 pu MW of un- 
contracted excess power in (11) i.e. 
&P L,uncoti= 0.05pu . Fig. 8 shows the stabilized 
transient response in the event of contract violation 
with the coordinated device operations. 
The convergence characteristic of the objective 
function given in (1) with the help of CRPSO algorithm 
is shown in Fig. 9. 

VI. Conclusion 

Eigenvalue analysis of an interconnected two area 
multiple unit hydro-hydro power system results in 
positive real parts for some of the eigenvalue pairs 
which confirms its dynamically unstable behavior. The 
optimized coordination of SMES-SMES and TCPS- 
SMES with the help of CRPSO effectively stabilize the 
system and significantly improves the transient 
response of area frequency and tie- line power 
exchanges under different contract variation 
inrestructured electricity market. The performance of 
SMES-SMES coordination yields the least value of 
Figure of Demerit and shows lesser undershoot and 
overshoot in dynamic responses of frequency and 
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Figure 9. Convergence characteristic obtained by CRPSO for TCPS- 
SMES and SMES-SMES coordination 

tie-line power exchanges as compared to TCPS-SMES 
coordination. 
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Appendix 

R=0.05,H=3s, D=l;a 12 =-l;T G =0.2s;T R =5s;T w = ls;R T =0.38; 

l+sT R 



B=2 l;Ti 7 =0.5 ;Gov=- 



l+sT G 



; Transient droop= 



l+s(R T /R)T R ' 



l-sT w 
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